Exactly solvable statistical model for two-way traffic 
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We generalize a recently introduced traffic model, where the statistical weights are associated 
with whole trajectories, to the case of two-way flow. An interaction between the two lanes is included 
which describes a slowing down when two cars meet. This leads to two coupled five-vertex models. 
It is shown that this problem can be solved by reducing it to two one-lane problems with modified 
C , parameters. In contrast to stochastic models, jamming appears only for very strong interaction 

between the lanes. 



I. INTRODUCTION 



'■^ The non-equilibrium properties of one-dimensional lattice gases have been studied intensively over the last years 

0. With lattice gases, one can model not only physical situations such as transport in solid ionic conductors Q or 
^ growth processes p, but also the traffic flow on roads Q|. Moreover, they can be used to study general features of 
I ' phase transitions in non-equilibrium systems For the traffic problem, the simplest model is the completely 

I asymmetric exclusion process (ASEP), where classical hard-core particles hop stochastically, with unit rate, in one 
direction only On a ring, one then finds a steady state of product form where all configurations are equally likely. 
In terms of the density p of particles, the flux is then given by j = p{l — p) and shows already the qualitative features 
d • found also in more sophisticated models, i.e. it vanishes for p = 0, 1 and has a maximum in between. 
^ An essentially new description of traffic flow was proposed recently by Brankov et al. In this work, non- 

' intersecting domain-wall lines on a square lattice were interpreted as space-time trajectories of cars. The weight of 
^ a trajectory is then obtained from the fugacities for horizontal and vertical moves. The single step, however, has 
Q . no stochastic interpretation. The problem can be formulated in terms of a five-vertex model which generates these 
lines and which is exactly solvable since it satisfies the so-called free-fermion condition. The result for the flux j is 
physically reasonable and very similar to that for a variant of the (stochastic) Nagel-Schreckenberg model |ll|] . In 
I ' this paper, we show that one can generalize this model to the case of two-way traffic where cars on different lanes 
^ interact with each other. The specific effect which we are treating is a tendency to slow down when another car 
is approaching. In the two-dimensional formulation, this is described by a modification of the fugacities whenever 
trajectories of oppositely-moving cars cross. One then is led to consider two five- vertex models with a certain coupling 
between them. It turns out, however, that this coupling only renormalizes the parameters in each subsystem, so that 
(«~>) . the problem remains solvable as before. One finds that, in this model, the effect of an obstacle, i.e. of a car in the 
' other lane, is relatively weak. While in stochastic models already a certain finite reduction of the hopping rate at one 
position usually leads to a traffic-jam phenomenon with a region of high density appearing in front of the bottleneck 
[0- jl^, this happens here only if the fugacity is reduced to zero for a large system. As will be explained, this feature 
is related to the different weighting of the trajectories in both cases. 

In the following, we first describe the model in section 2 and then explain its solution in section 3. Finally, in 
section 4, we discuss the results and add some further remarks. 



o 

O ■ II. MODEL 



We first recall the formulation of the original one-way traffic model in For a square lattice with periodic 

boundary conditions, the horizontal direction is interpreted as space, the vertical one as time (increasing downwards). 
Non-intersecting lines running towards the lower right are then drawn on the lattice and viewed as trajectories of 
right-moving cars. They do not end, so that the number A^i of cars is conserved. A horizontal step, representing a 
move, is given fugacity (weight) Xi, a vertical one fugacity ti. Statistical averages are then obtained from the partition 
function 
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(1) 



where Nrc{C) and Nt{C) are the total numbers of steps in the two directions for a certain configuration C of trajectories. 
These trajectories are generated with their correct weights if at each lattice site the vertices shown in Figure |l| are 
possible. 



FIG. 1. Vertex configurations for right-moving cars: Boltzmann weights are 0, l,xi,ti, \/x\t\ respectively 

Since crossings (vertex 1) are forbidden, one is effectively dealing with a five-vertex model which can be solved 
exactly via the Bethe ansatz, even for more general weights W5, wg p^ , p9[ |. In the present case, free-fermion 
techniques can be used to obtain the partition function ||20|] . 

For the two-way traffic model, we introduce a second lattice where trajectories run towards the lower left, corre- 
sponding to the cars in the other lane. The fugacities are taken to be t2 and the trajectories are now generated 
by the vertices in Figure 0. 



FIG. 2. Vertex configurations for left-moving cars: Boltzmann weights are 0, 1, X2,t2, \Jx2t2 respectively 



To formulate the interaction between cars in the two lanes, the indices at the vertices are specified in the following 
way 



lattice 1 



f3' 



lattice 2 



The variables a,/?, ... take the value 1 if a car is present (thick line) and zero otherwise. For all vertices, the so-called 
ice rule 



a + /3^a' + /3' and j + S = j' + S' 



(2) 



holds, which ensures the conservation law for the number of cars, separately for both lanes. 

We now imagine that the two lattices are placed above each other and attribute an additional Boltzmann weight 



(3) 



V = exp{—e) — exp ( {aS + a' 5' + [3"/' + [3'^) 



to adjacent vertices in the two layers. Then each crossing of two trajectories will be weighed with the factor 

< r exp{~h) < 1 (4) 
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To see this, one first notes that e = 0,v = 1 if one of the vertices is of type 2, i.e. if there is no car present. The 
values of e in the remaining cases are given in Table 1. It then follows that simple crossings, which involve a pair of 
vertices of type 3 and 4, lead directly to a factor r, see Table 1. 

Table 1. Interaction e between two adjacent vertices in the two layers. 
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If trajectories meet and run (anti)parallel before they separate again, each of the two branch points contributes a 
factor y/r. Some examples illustrating such crossings are shown in Figure |^. 
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FIG. 3. Trajectories of one left-moving and three right-moving particles with different types of crossings. 

One should mention that the choice (j^) for the interaction is not unique. The more general form for e 

e = A{a6' + p-f) + B{a5 + f3'-f) + C{a'5' -f + D{a' 5 + f3'i) (5) 

with A + B + C + D — h, still leads to the same factor r — exp{—h). The individual terms listed in Table 1, however, 
become more complicated. 

In the model defined in this way, one still has the freedom to choose the particle numbers and the fugacities. Thus, 
by setting X2 — 0, one can immobilize the cars in the second lane and treat in particular the case of one fixed obstacle, 
which is of special interest. 



III. SOLUTION 



We now show that the two-lane model can be solved by reducing it to the one-lane problem. The proof follows Ref. 
pl[ where a similar problem was treated. It is based on the ice rule (||) which relates horizontal and vertical bond 
variables. Suppose that the lattices have N columns and M rows, and let an,m{ln,m) and Pn,m{Sn,m) be the variables 
to the right and below the vertex {n,m), respectively, in the two layers. Then the total interaction is 

h *^ ^ 

m—1 n—1 
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With the help of (0), this can be rewritten as 



h I \ 
E=-^iYl ("0.™ + '^N,m)N2 + J2 ("^0,™ + 7JV,m)A^l + U (0) - U{M) (6) 

,m— 1 m—l / 



where 

N-l n 



U{m) = ^ ^(/3fe,m(5n+l,m - /?«+!, m5fe,m) (7) 



n=l fe=l 



contains only vertical bonds, while the other two terms in contain only horizontal bonds. Due to the periodic 
boundary conditions, the difference U (0) — U{M) vanishes and one obtains 



M M 



E = -hN2 ^ ajv,™ - hN^ ^ 7^,^ (8) 



m — 1 



This can be compared with the effect of a rescaling X\ xie X2 X2e which leads to an extra factor 
exp{—r]i{a + a')/2 — 772(7 + 7')/2) for each pair of adjacent vertices. Summed over all sites, this gives 



M N 
m—l n— 1 

which can be expressed as 

M M 

^' = -Yl Y^("0,m + aN,,n) + V{0) ~ V{M) - ^ y A^(70,r„ + lN,,n) " ■^(0) + V{M) 
m—l m—l 

where 

N-l n N N 

^(^") -EE -EE z^'^.- (9) 

n=l fc=l n=2 k=n 

and y(m) is defined analogously with f3 S. Using again the periodic boundary conditions, one finds 

M M 

E' = -rjiN aN,m - mN ^ 7JV,m (10) 

m — l m—l 

which has the same form as £' in (||). Therefore the interaction has the same effect as a change in the horizontal 
fugacities if one chooses ryi = hN2/N = hp2 and r]2 = hNi/N — hpi where pi and p2 are the densities of cars in the 
two lanes. The partition function is then 

Z{Ni,N2,Xi,X2,h,t2,r) = Z{Ni,xxrP\h)Z{N2,X2rP\t2) (11) 

This exact formula looks like the result of a mean-field treatment since only the densities in the other layer enter the 
expressions. One should point out that it also holds for more general choices of the vertex weights in the layers. Then, 
also the weight wi of vertex 1 has to be renormalized with the same exponential factor. 



IV. RESULTS AND DISCUSSION 



One can now make use of the results for the single-lane case 10 1. For one lane, the flux per site is equal to the 
average number of horizontal steps and given by 



NM 



1 



1 c 
—arccos — 
TT \ 1 



2x 



2xc - 



-1 -P 



(12) 
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where c = cos{ttp) and a; < 1 has been assumed. This is the physical region since the average speed of one car is 
V = x/{l — x). As a function of p, the flux has a maximum at p = {l/Tr)arccos{x), which shifts from p = 1/2 to p = 
as X increases. 

By inserting Xir^^ and X2r^^ into (14), one then obtains the fluxes ji and j2 in the two-lane case. These do not 
depend on the motion in the other lane, but only on the density there. Since j increases with x, the interaction factor 
rP" always reduces the flux, as expected. This reduction, however, becomes smaller as the density in the second lane 
decreases. For the case of only one car one has 

ji=j(pi,xiri/^) (13) 

and this approaches the value j{pi, xi) without interaction for large N. In order to slow down the traffic appreciably, 
one would need h N, i.e. an interaction increasing with the size, so that r vanishes exponentially. In other words, 
a transition occurs only at r = 0. 

As mentioned, the situation is different for stochastic models. There j shows a sudden decrease as soon as the 
corresponding quantity r (describing the reduced crossing probability at a defect) falls below a certain finite value re- 
This is connected with the appearance of a jam at the defect. In terms of trajectories, the effect can be described as 
follows. Consider a stochastic model as in |10(] where a particle can move an arbitrary distance horizontally, at each 
step continuing with probability p and stopping with probability q = 1 — p. At the defect, the quantities are p' < p 
and q' > q. A particle some steps away from the defect will typically move to the bottleneck and then stay there for 
some time. Due to q' > q, such a trajectory has a higher weight than any other one where it makes stops before and 
then crosses the defect immediately. The same holds for another particle following it, since this has q' = 1 once it has 
reached the site next to the first one. In this way, the jam builds as a region of vertical trajectories to the left of the 
defect. 

In the present model, the picture is different. There is no advantage in staying at the blockage, the crossing factor 
r and the weights x'^t^ are the same as for paths which approach the defect gradually. Nor is there an advantage for 
following particles to move next to the preceding one. Therefore no jam builds up. One could say that the model 
mimics the anticipation of disturbances by producing less densely packed trajectories. But, in shifting Vc to zero, it 
overestimates the effect. 

It is also interesting to compare the two models at the operator level. According to ||2^, the transfer matrix T of 
the (one-layer) five-vertex model commutes with the operator 

W = - E + «+i + 2^<) (14) 

n 

where H = (1 -I- a:^ — t'^)/2x, and it is easy to see that the ground state of Ti. gives the maximal eigenvalue of T. 
This operator shows very clearly the free-fermion character of the model and also its non-stochastic nature, since the 
necessary cr^cr^-terms (which are related to the loss processes in the master equation) are missing. 
If one uses more general vertex weights and wq, the operator 

W = - E (^n ^,t+l + (15) 



commutes with T, where A — (w^Wi — w^Wq) / {W2W4) |23|, p8|. Although this contains such terms and has the form 



of the time-evolution operator for fully asymmetric hopping |^ , , the fact that A is not equal to one still makes 
it different. On the other hand, this model is interesting, because it contains, in the x~t plane, a frozen phase with 
density p — 1/2 [Q, jl^, where the trajectories have the form of stairs with steps of unit lenght in both directions. 
This corresponds to synchronized traffic with always one empty site between the cars. As this phase gives the highest 
possible throughput of vehicles and persists for a wide range of parameters x, t, it represents the analogue of the 
maximal current phase in stochastic hopping models |26 ^. In the j-p relation, one then finds a cusp at p = 1/2. As 



mentioned above, also this model can be treated in the two-way case. However, apart from half-filling, the blocking 
properties will be similar to those described above. 
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